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ABSTRACT 


This report presents the results of calculations of the thermal radiation 
from the lunar surface incident onto a flat surface of unit area located a small 
distance above the moon. The orientation and height of the flat surface vary. 
The thermal radiation from the moon's surface is assumed to be non -diffuse. 

The calculations show that the thermal energy flux incident onto the flat 
surface can differ significantly for a lunar surface that emits radiation in a 
non -diffuse manner when compared to a lunar surface that emits in a diffuse 
manner. 
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CALCULATIONS OF NON-DIFFUSE INFRARED RADIATION 
FROM THE LUNAR SURFACE INCIDENT ONTO 
A UNIT ELEMENT ABOVE THE SURFACE 


SUMMARY 


The amount of infrared energy flux incident onto one side of a flat surface 
element of unit area located above the lunar surface has been computed. The 
orientation and height of the flat surface with respect to the moon's surface vary. 
The infrared emission from the lunar surface is not diffuse, but rather is 
emitted in a non -diffuse manner in accordance with a mathematical equation 
developed by Ashby. 

The results of the calculations show that for non -diffuse infrared emis- 
sion, contrary to diffuse emission, the orientation of the surface receiving the 
radiation and the sun position are important. 


INTRODUCTION 


Interest in the directional characteristics or non-diffuseness of the 
moon's surface infrared radiation has existed in this laboratory for several 
years; because of this, and to support the in-house studies, several contracts 
have been sponsored with outside research organizations. In one of these 
contracts, an empirical expression was developed by Ashby for predicting the 
infrared radiation emitted by the lunar surface. 

In the work presented here, a sunlit area on the surface of the moon 
has been assumed to emit non-diffuse infrared radiation, in accordance with 
the Ashby expression; and the amount of this radiation incident onto one 
side of a flat surface element of unit area located above the moon has been 
calculated. The orientation and height of the element with respect to the 
moon's surface vary. 



PURPOSE AND REVIEW OF WORK 


This research program in directional infrared radiation is an attempt to 
partially fulfill the goal of explaining and defining the lunar thermal environ- 
ment. One part of the directional radiation program is to develop a theoretical 
model with which one can accurately reproduce the existing experimental data 
and that is consistent with current knowledge of the physical structure of the 
moon's surface. Such a model should reveal important information about the 
physical features of the surface. Another part of the program is to present 
information about the directional aspects in a form useful to those solving 
problems in thermal experiments and thermal studies of the moon's surface. 

This report represents an effort to partially fulfill this last part. Calculations 
have been performed incorporating these directional effects; the results are 
graphically presented as infrared energy leaving the lunar surface and striking 
a unit area a small distance away. 

To support the in-house work, a number of contracts have been sponsored 
with outside research organizations. In 1965 a contract with Brown Engineering 
Co. led to an arrangement whereby one member of the Brown research staff 
spent 6 months at the Boeing Scientific Research Laboratories working closely 
with Saari and Shorthill. These two investigators have obtained an extensive 
amount of data over the past 5 years on the infrared emission from the moon 
[ 1] . Their measurements have been in the 10-12 p region of the spectrum. 

The purpose of the work between Brown and Boeing was to obtain from existing 
data a new set of measurements relating to the directional characteristics of 
lunar infrared emission. The effort consisted mainly of a reduction of many 
data points from 1 9 separate phases . The data were then presented graphically 
[2] in the form of brightness temperatures for various lunar locations and 
phases . A typical graph is shown in Figure 1 . 

Using the experimental results from the foregoing work and working 
under the sponsorship of Space Sciences Laboratory, Ashby and Burkhard have 
developed an expression for predicting the infrared emission that takes into 
account the directional behavior of the radiation [3] ; comparisons with experi- 
mental data show it to be reasonably accurate. The chief drawback of the 
formula is that it has an empirical rather than a theoretical derivation and, 
therefore, does not give any information about the physical characteristics of 
the lunar surface. 
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When the work was first begun, the only experimental infrared data that 
had been examined (to the author's knowledge) for any directional characteristics 
were those of Pettit and Nicholson [ 4] and Sinton [ 5] , and these data showed 
only a slight deviation from a Lambertian surface emission. It has been known 
for some time that sunlight is reflected by the moon's surface in a strongly 
directional manner. Although this is in a different part of the spectrum from 
the infrared, it did help to stimulate an initial interest in taking a closer look 
at the directional aspects of the infrared emission. The Pettit, Nicholson, and 
Sinton data and the photometric data were the basis for the speculation that the 
emission might be directional to a significant degree and deserved further 
investigation. 

The lunar landing of Surveyor I in 1 966 led to the conclusion that a 
pronounced directional emission did exist. The explanation of the temperature 
values recorded on the faces of the two instrument compartments of Surveyors 
I and in l 6, 7] required this directionality, and subsequent investigations of 
earth -based measurements by Saari and Shorthill [8] reaffirmed this. 

Some of the Saari and Shorthill data showing the directional aspects are 
presented in Figure 2. Measured results from earth -based telescopes are 
compared with calculated results using the previously mentioned Ashby expres- 
sion. For each elevation angle of the sun, the brightness temperature is shown 
(for a diffuse surface the brightness temperature is a single point and is shown 
as an open circle) . By comparing the diffuse value with the non -diffuse points, 
the degree of directionality can be ascertained. Notice that the radiation exhibits 
the greatest non -diffuseness, i.e. , directionality for small sun elevation angles. 
Another point to notice is that the largest temperature values occur, for each 
sun elevation angle, when the surface is viewed from the sun direction: that is 
when an observer views the surface with the sun behind him . This means that 
a greater part of the infrared radiation is emitted toward the source, i.e. , the 
sun, than in any other direction. Such behavior may be referred to as back- 
emitted, a term analogous to backscatter for the reflected sunlight from the 
moon that is reflected or scattered more toward the source than in any other 
direction. 

Gross surface roughness has been suggested as the reason for the back- 
emitted directional radiation. This thought was put forth by Pettit and Nicholson 
to explain their measurements of non -diffuse radiation, and still appears to be 
the most satisfactory explanation. The surface, if assumed to consist of peaks 
and valleys (or of numerous craters and rocks) , will receive sunlight mainly on 
the side facing the sun. The other side will be in the shadow and will, therefore, 
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be at a much lower temperature . This condition will exist to the greatest 
degree when the sun is near the horizon, and to the least degree (or not at all) 
when the sun is straight overhead. 

The earth -based measurements were made with telescopes that have 
resolutions of the order of 8 - 10 seconds of arc. At the lunar distance, this 
means that the measurements average the radiation (at the moon's center) 
over a circular region of approximately 14-18 km in diameter. Whether the 
measurements would reveal this same directional characteristic over a much 
smaller region, say on the order of a few meters, is an unanswered question. 
The surface would probably be more diffuse. The reason given earlier — gross 
surface roughness — would not suffice to explain any directional behavior 
•observed on this scale. 


THERMAL EMISSION PER UNIT SOLID ANGLE 


The emission of infrared energy from a local area on the moon in a 
given direction can be computed from Lambert's equation, I( e) = I 0 cos e, if 
the emission is assumed to be diffuse. However, if the emission is non-diffuse, 
as recent data indicate [7] , the infrared emission from a local area on the 
moon in a given direction can be more accurately determined from the following 
expression, which was developed under a NASA contract [3] . 


T/ . _ , a^osi + ajcosa* a, . , . , v 

I 0 (i, e, a) = [ (7T - f of |) cos | a | + sin ) or | ] 


1 + a. 


sin a 
cos i 


( 1 ) 


where: 


a = — 


i 2 + e 2 - 2ie cos( - <f > ) 

2 ^ 7r 2 4i 2 e 2 

J + ~ 2ie eos( 0. - <p ) 


I 0 (i, e, a ) = Infrared radiance in watts/m 2 - steradian 
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at = 335.417 
a 2 = 97.626 
a3 = 51 . 603 
a 4 = 84.377 

(All of the a's are empirically determined using the Saari 
and Shorthill data of references 1 and 2, and have units of 
watts/m 2 - steradian. ) 


i ,e, 4> <p a = Angles defined in Figure 3. 

This expression holds true only for a sunlit surface; therefore, it is not useful 
during the lunar night or for lunar areas that are completely in shadow. From 
equation (1) the infrared energy is 

I(i, e, a) = I 0 ( i, e, a) cos e (2) 

The brightness temperature of a local lunar area when observed from 
different directions and for different sun elevation angles can be obtained from 
equation ( 1) as follows: 

I 

^ h a ) 4 

a (3) 

where: 



a = Stefan -Boltzmann constant (5.673 x 10 -8 watts/m 2 - °K 4 ) 

The brightness temperatures as obtained from equation ( 3) for two sun elevation 
angles are shown in Figure 4. 


TOTAL THERMAL EMISSION 


If the infrared emission in the direction of dA 2 from many local areas, 
such as dAj (Fig. 3) , is summed over a local lunar surface area A t , the energy 
incident onto dA 2 can be determined. The energy from dAj to dA 2 is found by 
using equation ( 4) below. The angles e and d are required to compensate for 
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the skewness of the surfaces with respect to each other, while the term r 2 repre- 
sents the attenuation of the energy in accordance with the inverse square law. 


dE 

where: 


In( i, e, a;) cos e cos 0 
— 

r * 1 


dAj dA 2 


(4) 


dE is the energy flux from dA t to dA 2 in watts 

The lunar area A t is taken as flat and circular with a radius R. Summing 
over A t and assuming dA 2 to be unity gives 


E _ r I n (i, e, a) cose cose ^ (5) 

Ai 

where: 

E is the energy flux from A t to dA 2 in watts/m 2 
Formulating r and dA t in terms of the fundamental angles gives 

dAj = h 2 tan e sec 2 e d $ d e (6) 

r = h sec e ( 7) 

Substituting equations (6) and (7) back into equation (5) gives 
tan -1 r/H 27 r 

E = f f I 0 (i, e, a) cos 0 sin e d <f> d e (8) 

J o 

The angle 6 can be obtained for any orientation of dA 2 and in terms of e, y, 
and <p by 


9 = cos -1 [ cos e cosy - sin e cos( 180 - <f>) siny], 0 < 7r/2 

( 9 ) 

9 = 7t/2 , 0 > 7r/2 
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Because of the limits of integration, the integration over <p is done in four 
parts — one part for each quadrant of A t . Finally, 


tan -1 R/h 0** 

E = f f I 0 ( i, e , a) cos 0 sin e d <p d e 

‘'o •'o 

tan -1 R/h v - <t>* 

+ f f I 0 ( i , e , a ) cos 0 sin e d <p d e 

0 tt/2 

tan -1 R/h 3ir/2 

+ f J I 0 (.i, €, Qf) cos 0 sine d $ d e 

0 7T + <f>* 

tan -1 R/h 27r 

+ / / lo (i» €, a) cos 0 sin e d <p d e 

0 37r/2 + <£*** 


( 10 ) 

The values for , <p'' :< * , and * * depend on the orientation of dA 2 . The 
values and corresponding orientations are as follows: 


0 

<P >] 

<1 P* 


* 


= 0 


= cos 


-if tan |y - 90 | 
tan e 


= 7r/2 


= 7r/2 

(j)** = o 


j J tan | y — 90 1 


= cos' 


tan e 


ou 

■V* ^ ^ 




^ 

# ^ 

<P 


= o 


<$> 


>',< >I< 5 ( '< 


= sin 


= tt/2 


-l / tan j y - 90 | 


tan e 


y<ir/2, 0 < e — j y - 90 j 
y < 7r/2 , | y - 90 | < e < tan -1 R/h 
y -- tt/2 

y ^ 7r/2 

y > 7t/2 , |y - 90 | < e 
y > 7r/2 , |y - 90 1 s e 

y ^ 7r/2 

y > tt/2 , I y - 90 1 < e < tan -1 R/h 

y > V 2 
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RESULTS OF COMPUTATIONS 


Equation (10) has been integrated numerically. Figure 3 provides an 
explanation of the angles involved in the integration, and Figure 5 is helpful in 
understanding the results. Figures 6-17 present the results of the computa- 
tions using equation ( 10) . Shown on each figure is the infrared energy flux 
from an area Aj of the lunar surface (assumed to have an infrared emittance 
of unity) incident onto one side only of the element A 2 for different sun eleva- 
tions. Each curve represents a different tilt angle y for A 2 . Each group of 
figures (four to a group) represents a single azimuthal orientation 0 2 , of 
the element. When <p 2 is 180°, the element faces west. These curves are 
intended to have usefulness primarily for purposes of thermal design. 

It is instructive to delineate as plainly as possible the directional aspects. 
One way of doing this is by comparing a diffuse and non-diffuse case for one 
orientation of the element, as is done in Figure 13. As in the previous figures, 
energy is plotted against sun position. A better way of showing the directional 
characteristics is by normalizing the energy values in Figures 6-12 about the 
maximum value (the curve for R/h = 1000, y = 0) for each azimuthal 
orientation, <p 2 , of the element. This normalized result or view factor is 
more useful for performing thermal radiation calculations. Figures 14-20 
show the view factor presented as in the previous graphs. The difference in 
the figures, as before for the energy curves, is the azimuthal orientation of 
the element and the ratio of size of the lunar surface to the height above the 
surface of the element. Figure 14 is typical of the other figures and will be 
used to point out some of the significant points. In Figure 14 the element 
always points west, i. e. , <p 2 = 180° . As the sun moves from the sunrise to 
sunset position, the directional aspects of the radiation can be plainly seen on 
comparison with the diffuse value, which is shown as a single point (although 
it could be shown as a horizontal line through this point) . Notice particularly 
the general shape of the curves, how the energy striking the element is greatest 
near sunrise, i.e. , when the sun elevation angle is small and decreases as the 
sun moves from this position through lunar noon, finally having the smallest 
view factor when the sun is near sunset and facing the element. For an element 
facing east ( cp 2 = 0°) , the shape of the curve is exactly reversed, i. e. , the 
smallest view factor occurs at sunrise and the greatest occurs at sunset. The 
other figures are for elements facing in other directions. Notice the element 
facing north (<p 2 - 90) , Figure 20. The curves are almost flat, meaning that 
almost no changes in the surface radiation occur for changing sun angle when 
the lunar surface is observed from such a direction. This is somewhat evident 
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from Figure 4 where the shape of the model indicates the distribution of the 
energy. Notice that the radiation is asymmetrically in the plane containing 
the surface normal and the east -west axis (the plane containing the sun vector) 
and symmetrical in the plane containing the surface normal and the north -south 
axis. 


The curves shown in Figures 14-20 are for an element confined to one 
quadrant of the azimuthal plane, i.e. , <p 2 = 90° - 180°. For elements facing 
directions other than this, symmetrical considerations depicted in Figure 4 can 
be used to get the view factor from the values shown here. For example, the 
view factor values are the same for <f> 2 = 225° as for <p 2 = 135° . Figure 4 
aids in understanding why this is so since the figure shows that the brightness 
temperature or energy in these two directions is the same. For <p 2 = 45° , the 
values are the same as for <p 2 = 135° , but the abscissa on the graph must be 
reversed, i.e., the scale should begin with 180° and end with 0° . When this is 
done, the view factor values that were at 30° sun elevation (for example) will 
appear instead at 150® sun elevation. This reversing scheme must always be used 
when the element faces between = 0* - 90* or 270* - 360* . 

The curves clearly show the effects of the back -emitted radiation since 
the largest view factor occurs when the element faces away from the sun. 

This variation of energy with element position can be shown in another 
way by having the element rotate about an axis vertical to the moon's surface 
(i.e. , rotation in the azimuthal plane) while the sun's position is fixed. This is 
shown in Figure 21. The element itself is vertical to the surface (y = 90°) 
and the diffuse values are again shown as a single point. Notice that the radia- 
tion shows no change when the sun is straight overhead ( c o - 90° ) . This is 
because the radiation is symmetrical with respect to the azimuthal plane. This 
trend toward symmetry in the azimuthal plane as the sun approaches lunar noon 
can be observed in Figures 2 and 4. Infrared measurements of the full moon, 
such as those of Pettit and Nicholson, are for this condition; and directional 
aspects when viewed from earth are at a minimum. 


FINAL REMARKS 


One method for taking the directional aspects of the lunar surface infra- 
red radiation into account in numerical calculations has been shown. The 
results presented here have shown that calculations based on diffuse radiation 
may be as much as 80 W/m 2 ( 16 percent of total) at variance with calculations 
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based on non-diffuse radiation (Fig. 13) . Orientation of the surface receiving 
the radiation and sun position are also important. For a sun elevation angle of 
30°, the radiation striking a surface which is oriented vertically to the lunar 
surface may vary as much as 120 W/m 2 (33 percent of total) as the surface is 
rotated 180° about a vertical axis (Fig. 21) . 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, December 22, 1967 
129-04-02-00-62 
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(The elevation angles are measured from the lunar surface 
in the direction of sunrise and in the plane containing 
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(a) Angles Involved in Integration of Equation ( 10) 



(b) Angles Used in Equation (1) 
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FIGURE 3. ANGLES INVOLVED IN DISCUSSION AND CALCULATIONS 
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indicates the magnitude of the temperatures) 





SUN ELEUATION-OEOREES SUN ELEVATION-DEGftEES 


(a) 0 2 = 180°, R/h = 1000 


(b) 0 2 = 180°, R/h = 10 




(c) 02 = 180°, R/h = 2 (d) 0 2 =18O O , R/h = 1 


FIGURE 6 . ENERGY INCIDENT ONTO ONE SIDE OF UNIT AREA (dA 2 ) 
VERSUS SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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ENERGY (wotU/m*) ENERGY (wom/m*) 



(c) 0 2 =15O O , R/K = 2 


(d) 0 2 = 150° , R/E= 1 


FIGURE 8. ENERGY INCIDENT ONTO ONE SIDE OF UNIT AREA (dA 2 ) 
VERSUS SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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SUN ELEVATION -DEGREES 


(d) $ 2 = 135°, R/h = 1 


ONE SIDE OF UNIT AREA (dA 2 ) 
VARIOUS ORIENTATIONS OF dA 2 
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(c) 0 2 = 120°, R/h = 2 


(d) 0 2 = 120° , R/h = 1 


FIGURE 10. ENERGY INCIDENT ONTO ONE SIDE OF UNIT AREA (dA 2 ) 
VERSUS SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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FIGURE 11. ENERGY INCIDENT ONTO ONE SIDE OF UNIT AREA (dA 2 ) 
VERSUS SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 
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FIGURE 12. ENERGY INCIDENT ONTO ONE SIDE OF UNIT AREA (dA 2 ) 
VERSUS SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA, 
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FIGURE 14. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA z ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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(d) <^2 = 165°, R/h = 1 


FIGURE 15. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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FIGURE 16. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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FIGURE 17. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 












FIGURE 18. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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FIGURE 19. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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FIGURE 20. NORMALIZED ENERGY (VIEW FACTOR) FOR RADIATION 
BETWEEN LUNAR SURFACE AND ONE SIDE OF UNIT AREA (dA 2 ) VERSUS 
SUN ELEVATION ANGLE FOR VARIOUS ORIENTATIONS OF dA 2 
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